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We investigate the influence of cosmic magnetic fields on gravitational wave perturbations, and find 
exact solutions on large scales. We show that a large-scale magnetic field can generate large-scale 
non-decaying gravitational waves. In the general case where gravitational waves are generated by 
other mechanisms, a large-scale magnetic field introduces a new decaying tensor mode and modifies 
the non-decaying mode. The direct effect of the magnetic field is to damp the gravitational waves, 
while an indirect magneto-curvature effect can either damp or boost the waves. A magnetic field 
also leads to a breaking of statistical isotropy, and the magnetic imprint on the tensor spectrum in 
principle provides a means of detecting a primordial field. 
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I. INTRODUCTION 

Magnetic fields seem to be everywhere that we can look 
in the universe, from our own sun out to high-redshift 
Lyman-a systems. The fields we observe (based on syn- 
chrotron radiation and Faraday rotation) in galaxies and 
clusters have been amplified by gravitational collapse and 
possibly also by dynamo mechanisms. They are either 
primordial, i.e. originating in the early universe and al- 
ready present at the onset of structure formation, or they 
are protogalactic, i.e. generated by battery mechanisms 
during the initial stages of structure formation. One way 
to distinguish these possibilities would be to detect or 
rule out the presence of fields coherent on cosmological 
scales during recombination via their imprint on the cos- 
mic microwave background (CMB) radiation. Consider- 
able work has been done to investigate the nature of the 
magnetic imprint on the CMB (see jjj for recent reviews 
and further references). 

We show here that cosmological magnetic fields also 
leave a characteristic imprint on the cosmological gravi- 
tational wave background. Direct detection of this back- 
ground, especially on large scales, is not likely for a 
considerable time (see, e.g., ||), so that our results do 
not provide a practical means for detecting or limiting 
a large-scale magnetic field. However, these results are 
a necessary first step in a theoretical understanding that 
can be developed in anticipation of gravitational wave de- 
tection. Furthermore, our results suggest a similar inves- 
tigation of the magnetic imprint on gravitational waves 
generated by compact objects. This astrophysical prob- 
lem is more complex than the cosmological problem that 
we discuss, but it is likely to lead to larger effects with 
stronger prospects for observational detection. 

The effects of gravitational waves on electromagnetic 
fields, in particular, the generation of electromagnetic 
pulses by gravitational waves, has been previously con- 
sidered, including the implications for gravitational wave 



detection || (see Q{ for recent results and further refer- 
ences). Here we consider a different and new effect, i.e. 
the influence of a large-scale magnetic field on gravita- 
tional waves, including the capacity of the field to gener- 
ate gravitational waves. This latter possibility has been 
previously investigated on small scales in || , where tan- 
gled magnetic fields during recombination were shown to 
generate gravitational wave perturbations, which in turn 
contribute to CMB anisotropies, thus providing a way of 
limiting the strength of the magnetic field. 

We consider in a general cosmological context the prob- 
lem of how magnetic fields affect gravitational waves, fo- 
cusing on large-scale fields. A large-scale homogeneous 
magnetic field is strongly limited by large-angle CMB 
anisotropies |J: B now < 1 0~ 9 G. Although the maxi- 
mal energy density supported by the field is a very small 
fraction of the total cosmic energy density, this is also 
the case for gravitational waves, so that in principle, as 
we show, the magnetic effect on these waves need not 
be negligible. We find exact analytic solutions on large 
scales. These solutions show that the direct effect of the 
magnetic field is to damp the gravitational waves, while 
an indirect effect due to the coupling between the field 
and the spatial curvature can either damp or boost the 
waves. Qualitatively, the presence of a large-scale cosmic 
magnetic field is signalled by a breaking of global statisti- 
cal isotropy in the gravitational wave spectrum, which is 
in principle detectable (anisotropy detection is discussed 
in general in H]). 

II. PERTURBATION EQUATIONS 

The background is a spatially flat Friedmann space- 
time with a non-magnetized perfect fluid, representing 
non-interacting baryons, radiation and cold dark matter 
with the same 4-velocity u a , and Hubble rate H = a/a. 
The background equations are 
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-3H(p + p) 



3H = p + A, 
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where p and p are the total energy density and pres- 
sure. The perturbed universe is permeated by a weak 
large-scale magnetic field B a , whose energy density 
p mag — \B 2 is a first-order quantity (so that B a is 
"half-order" ||). The magnetic field also has first-order 
isotropic pressure \B 2 and tracefree anisotropic stress 

TTah = —B( a B b ) = -[/i( a c /ib) d — \h cd h a b\B c B d , (3) 

where h ab = g ab + u a u b projects into the comoving rest 
space and the round brackets on indices denote sym- 
metrization. The baryonic fluid is treated as a perfect 
fluid of infinite conductivity, with energy density ph and 
isotropic pressure pb- High conductivity ensures that 
any electric fields that might have been present dissi- 
pate quickly. For tensor perturbations, which do not ex- 
cite relative velocity perturbations, u a is the 4- velocity of 
baryons, radiation (p r = \pi) and collisionless cold dark 
matter. 

Both the isotropic and anisotropic magnetic effects are 
of the same perturbative order as gravitational waves. 
The covariant Maxwell equations Q_ imply the induction 
equation, whose nonlinear form is 



h a B b — (a a i, + e a bcU c — %Oh a b) B , 



(4) 



where cr a b, to a and O are respectively the shear, vorticity 
and expansion. In the background, = 3H and uj a — 
= (Tab- The induction equation implies the linearized 
conservation equations P,|10[ 



(B 2 )' = —AHB 2 



where the dot denotes it a V Q . It follows that 
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TTafc = -B 2 (— ) n {a n b) , n a = , n a = . (7) 

The unit magnetic direction vector n a is parallel- 
propagated along each observer world-line to first order. 
Maxwell's equations also imply the constraint [pTj 



D 6 7r a6 = e abc B b cui\B c - ±D a B 2 



(8) 



where D a is the projected covariant derivative: 
DaS "c-.- = (V a S b '" C ...)_L, which leads to covariant spa- 
tial curl operators: curlS^ = £ a fccD S c , curl S a b = 



~cd(a 



c Qd 



D C S" 
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The 4-acceleration A a is determined by Euler's equa- 
tion |@ 



(Pb +Pb)A a = -DaPb + e abc B b c\xA B c . 



(9) 



For tensor perturbations, the 4-acceleration should van- 
ish at first order. The non-magnetized condition for 
transverse traceless (pure tensor) modes is ]l3| , |l2t 



= D fI p = D a p. 



(10) 



When there is no magnetic field, this ensures that A a = 
= D a 9 and that all the traceless tensors, such as a ab , 
are transverse. In the magnetized case, two additional 
constraints need to be imposed: 



D a B 2 = = e abc B cmlB c 



(11) 



i.e., the magnetic energy density is homogeneous and the 
field is force-free, to first order. Equations ( |Io| ) and ( |Tl] ) 
together imply 



A a = = D„0 , 



(12) 



and that all traceless tensors are transverse. (The ex- 
pansion gradient is seen to vanish from the propagation 
equation for D a/ o S.) The magnetic field is felt only 
through its transverse traceless anisotropic stress n ab . 

Gravitational radiation is covariantly described by the 
electric (E ab = £7( a &>) and magnetic (H ab = i?( a &>) parts 
of the Weyl tensor, which support different polarization 
states fl^| , and which obey equations remarkably analo- 
gous to Maxwell's IIJQ: 



E ab = -3HE ab + ^Hn ab - \ (p + p)a ab + curl H ab , (13) 

(14) 
(15) 
(16) 



H ci b = -3HH ab - cmlE ab + icurl7r ab , 



0~ab — —2Ha a b 
Hab = curl a ab , 



where 



V b E ab 



B b H n 



D b a ab 



D b 7r ab 



(17) 

In the magnetized case, ir a b is given by Eq. (Q). 

The direct relation between the shear and the mag- 
netic Weyl tensor means that we can describe gravity 
wave evolution via E ab and a ab alone. Equations (|l^)- 
( |l5| ) show how the magnetic field is a source of shear and 
gravitational wave perturbations via its anisotropic pres- 
sure. Equation ( |l6| ) allows us to decouple Eq. ([l4]) from 
the system, which reduces to Eq. (|l5| ) and 
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3HE ab + iHitab -Up + V)Oab - D 2 CT ab , (18) 



on using the linearized identity curl curl S a b — —^Sab- 
It follows via the role of 7r ab in these equations that there 
is a directional effect at all scales on gravitational waves 
due to the magnetic field. We will further discuss this 
effect on large scales below. Note also that to first order, 
the evolution of the magnetic field is independent of any 
gravity wave effects. The field affects gravity waves, but 
there is no backreaction on the field. 

We can further reduce the system of equations to a 
single covariant wave equation in a a b'- 



V 2 a ab - (T a b = 5H&ab + [| (p- 3p) + 2A] a a b 



+ 2B; 



(19) 
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on using Eqs. (g), (§), (@), (||) and ©. The solution of 
this wave equation then gives the gravito-electromagnetic 
tensors via Eqs. ( filf ) and ([l6|). In practice, it is better to 
solve the first order coupled system for E a b and o a b than 
the single second order wave equation. 

The covariant Maxwell- Weyl approach to gravitational 
wave perturbations may be related to the metric-based 
approach (see also |jlj|[l^,[l4|). The transverse traceless 
metric perturbation is defined by 



ds 2 



-dt 2 + a 2 [S, tj + f ij ]dx i dx j 



(20) 



In these coordinates, it follows that E na = and |L5[ 



E; 



fij + Hfa 



(21) 



on large scales. The wave equation for /y then shows 
that, neglecting anisotropic stresses, = \Hfij. Thus 
the energy density in gravitational waves, p gw — \fijf 1 ^, 
is given in the absence of a magnetic field in Maxwell- 
Weyl form as 



Pgw 
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T7i Tpab 



(22) 



on large scales. Whether or not there is a magnetic field, 
the dimensionless contribution to the power spectrum on 
large scales is determined by X 4 E a i,E ab , where A is the 
wavelength [fl2f . Thus a dimensionless measure of ampli- 
tude on large scales is given by 

T 2 = a A E ab E ab . (23) 

III. LARGE-SCALE SOLUTIONS 

On superhorizon scales we can neglect the Laplacian 
term in Eq. (|l8|). In the radiation era, Eqs. ( |l5| ) and ( |l8| ) 
then have the solution 

E ab = C ab T -l + C ab T -5/2 _ l^b) ^-2 ^ ^ 

H a ab = 2Cf - C ,ab T- 3 ^ - |n< n 6 >Bgr -1 , (25) 
where CJ. = an d r = t/to. In the matter era: 

E ab = C ab T -4/3 + ^-3 _ 3^)^-8/3 ^ (2g) 

H a a ab = C ab r- X ^ - §C™ b r- 2 - in< Q n b >S 2 r- 5 / 3 . (27) 



The magnetic field introduces a new mode of gravita- 
tional wave perturbations on large scales, which decays 
less rapidly than the standard, non-magnetized decaying 
mode, and it modifies the standard non-decaying mode 
(of a 2 E a b), i-e., magnetic terms enter C" b (see below). 
The directional influence of the large-scale magnetic field 
on gravitational waves means that the field breaks the 
statistical isotropy of the large-scale tensor spectrum, as 
further discussed below. 



The new magnetized modes mean that if there are no 
gravitational wave perturbations initially present, then 
these can be generated by a large-scale magnetic field. 
(The generation of tensor perturbations by small-scale 
tangled magnetic fields is investigated in Ref. H.) This 
could happen if a large-scale magnetic field is created at 
time to. Large-scale magnetogenesis can occur for ex- 
ample in the recombination era [ fl6[ or the preheating 
era after inflation (see Jl7j] and references cited therein). 
Since E ab (t Q ) = = E ab (t ), it follows from Eq. (§|) that 



E a b — 



-5/2 



±T" 2 
2 ' 



B 2 n 



(a«b> 



(28) 



if to is in the radiation era, with a similar result if to 
is in a matter-dominated era. These 'purely magnetic' 
gravitational waves have a non-decaying mode, 
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r - -Jg - a o B o ■ 



(29) 



and they satisfy E a b oc n ab , so that in some sense they 
are maximally anisotropic (see below). The generation 
of tensor perturbations via magnetogenesis may be com- 
pared with the magnetic generation of density perturba- 
tions, which produces a growing mode Jl8| , ^0[ . 

In the general case, tensor perturbations are generated 
by other mechanisms and then influenced by the mag- 
netic field. The dimensionless measure of amplitude, T, 
defined in Eq. (f||, is seen from Eqs. to be 

constant at late times (t ^> to), i.e., when we neglect the 
decaying modes. Explicitly, we find the following mag- 
netic influence on the non-decaying mode of large-scale 
gravitational waves. 
In the radiation era, 

r 2 = I [4r 2 + 4a 2 H 2 Z 2 - 8a Q H Q Xo] 

+ §f a l B2 o ( a l B o - 6£ o + 6a H S ) , (30) 

where we expressed C? in terms of the dimensionless 
physical scalars, 

S 2 = a 2 Oab° ah , S = aaabn a n b , (31) 
£ = a 2 E ab n a n b , X = a 3 E ab (J ab . (32) 

In the matter era, 

r 2 = — [r§ + a^Ho^l - 2a HoXo] 

+ ^-a 2 B 2 (a 2 B 2 - 6£ + 6a H So) . (33) 
50 

The effect of the magnetic field is more clearly brought 
out if we use the Gauss-Codazzi equation || 

B (ab) = B ab — Ha a b + \^ab , (34) 

where R* ah is the Ricci tensor of the spatial hypersurfaces 
orthogonal to u a . This equation implies 
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K = £- aHS 



\a 2 B 2 



where 



K = a 2 R* {ab) n a n b 



(35) 



(36) 



is a dimensionless curvature scalar, giving the anisotropic 
3-Ricci curvature along the magnetic field, due to the 
gravitational waves. Using Eq. (|35|), we can rewrite 
Eqs. ph and (|3T 



r- 



1 



4a 2 l _ff ( 2 £ 2 l — 8a HoXo] 



2 
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— a Q B Q ( a B 



in the radiation era, and 



r 2 = I \r 2 



25 



6^ ) , 



(37) 



(38) 



in the matter era. 

These equations show that there are two aspects to 
the magnetic effect: a 'pure' magnetic effect (neglecting 
curvature), and a magneto-curvature effect. The pure 
effect, proportional to — Bq, serves to damp the ampli- 
tude relative to the non-magnetized case. It is due to the 
tension of the magnetic field lines, which means that the 
magnetic field resists the distortions induced by a grav- 
itational wave yj. The magneto-curvature effect, pro- 
portional to —BfilZo, arises from the geometric coupling 
between curvature and the field in general relativity, due 
to the vector nature of the field |p| [TT] , p^| . It depends 
on the sign of TZq, i.e., on the sign of the anisotropic 3- 
Ricci curvature component along the magnetic direction 
at time to. If IZo > 0, then the damping is reinforced; 
otherwise, the magneto-curvature effect enhances the am- 
plitude, acting opposite to the pure magnetic effect. 



IV. MAGNETIC IMPRINT ON GRAVITATIONAL 
WAVES 

In the above solutions, the terms in round brackets give 
the magnetic correction to the non-magnetized model. A 
rough estimate of the relative magnitude of the magnetic 
imprint on the large-scale gravitational wave spectrum is 
therefore 



r 2 



r 2 | B =o 



a B 
1 



(39) 



Now Pmag/Pr is constant to first order, and large-angle 
CMB anisotropies at the 10 -5 level imply the rough limit 
(Pmag/Pr)doc < 10~ 5 . The CMB quadrupole also places 
an upper limit on large-scale tensor perturbations [B0[ : 
(O gw /n r ) now < 10 -10 . Clearly a is sensitive to the ten- 
sor contribution to CMB anisotropies: the smaller this 
contribution, the higher a is. However, this increase 
in a comes at the cost of greater difficulty in detect- 
ing the gravitational wave background. In the most 
favourable case for detection, when the tensor contribu- 
tion to the CMB quadrupole is of the same order as the 
scalar contribution and the magnetic contribution to the 
CMB is maximal, then, with to — idec, it follows that 
a < (p r /p)dec ~ 10 _1 . 

Thus in principle, observations of the large-scale gravi- 
tational wave background can provide quantitative limits 
on a possible large-scale magnetic field, although in prac- 
tice, detection of such a background presents great dif- 
ficulties (see and references therein). Thus the CMB 
remains the best way to limit a magnetic field, on all 
scales |||||||. 

Qualitatively, the presence of a magnetic field is sig- 
nalled by the breaking of global statistical isotropy in the 
tensor spectrum. This may be compared with a simi- 
lar situation regarding the effect of a large-scale mag- 
netic field on CMB polarization anisotropy, where a mag- 
netic presence is signalled by a correlation between B- 
polarization and temperature anisotropies, which is im- 
possible if statistical isotropy holds . A possible mea- 
sure of global statistical anisotropy is the average over all 
directions n a of the quantity 



r 2 - 3£ 2 

f 2 



(41) 



In the absence of a magnetic field, statistical isotropy will 
give (£) n = 0, since 



fc\2\ 



((E ab n a n») 



ab 



E ab 



(42) 



A measurable magnetic field will produce a non-zero £ for 
the magnetic direction n a . In the case of gra vitational 
waves generated by a magnetic field, Eqs. (|2|) and (|4l| ) 
show that |£| = 1. 

Although we were only able to find exact solutions on 
large scales, the small-scale solutions should share the 
same basic features. Finally, we point out that if the same 
qualitative physical effects of magnetic fields on gravita- 
tional waves occur in the context of compact objects, 
then since the gravitational and magnetic fields can be 
very strong, this could have more immediate and impor- 
tant observational implications. This is a topic currently 
under investigation. 



On using Eqs. 



P 2 

rmag 



23) and (E3J), we find 



(Pmag/ Pr)p 
(Pgw/Pr)o 



(40) 
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